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We  demonstrate  the  effects  of  shear  modulation  on  the  viscoelastic  response  of  nanometer-scale  single-asperity  contacts  under  static 
and  dynamic  loading  conditions.  Contact  stiffness  and  relaxation  time  are  determined  for  contacts  to  poly(vinylethylene)  using  a  scanning  force 
microscope  (SFM).  Knowledge  of  the  torsional  stiffness  K©  of  the  SFM  cantilever  is  not  required  to  determine  the  relaxation  time.  The  relaxation 
time  was  several  orders  of  magnitude  slower  than  the  bulk  relaxation  time  but  decreased  slowly  to  the  bulk  value  as  the  sample  age  increased. 
Contacts  showed  no  evidence  of  microslip.  We  show  that  the  shear  response  observed  during  the  making  and  breaking  of  the  contacts  provides 
information  about  the  time  evolution  of  the  contact  area  that  is  not  available  in  force  vs.  distance  curve  measurements. 
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1.  Introduction 

Nanometer-scale  single  asperity  contacts  have  broad 
importance  in  tribology.  They  are  a  key  element  in  micro¬ 
scopic  models  of  macroscopic  contacts  as  formulated  by 
Bowden  and  Tabor  [1],  their  sizes  are  comparable  to  me¬ 
chanical  contacts  in  MEMS  devices  [2],  and  they  are 
widely  used  as  models  of  AFM  contacts  [3].  In  each  of 
these  areas,  it  is  essential  to  understand  the  mechanical 
behavior  of  the  contact,  particularly  under  combined  nor¬ 
mal  and  shear  loading.  In  this  letter  we  describe  new  ex¬ 
periments  using  the  scanning  force  microscope  (SFM)  to 
demonstrate  the  effects  of  viscoelasticity  on  the  mechanical 
properties  of  single  asperity  contacts.  The  possibility  of 
making  these  measurements  was  suggested  previously  [4], 
The  method  is  based  on  well-developed  mechanical  modu¬ 
lation  techniques  used  to  measure  the  viscoelastic  mechani¬ 
cal  properties  of  bulk  polymers  [5]  and  has  been  used  pre¬ 
viously  to  study  elastic  behavior  in  SFM  contacts  [6-8]. 
Similar  methods  have  also  been  applied  to  much  larger 
diameter  contacts  using  the  surface  forces  apparatus  [9,10] 
and  microscopic  sphere-on-flat  contacts  [11].  The  tech¬ 
nique  is  also  closely  related  to  friction  loop  methods  [12] 
and  to  SFM  methods  that  use  modulation  normal  to  the 
sample  surface  [13,14]. 

Figure  la  illustrates  schematically  how  the  shear  re¬ 
sponse  is  measured.  An  SFM  tip  of  height  H,  mounted  on  a 
cantilever  force  sensor,  is  initially  in  contact  with  a  vis¬ 
coelastic  sample  above  O  under  applied  load  FN.  The  sub¬ 
strate  is  displaced  parallel  to  its  surface  by  distance  xa  to 
O’.  This  causes  a  displacement  of  the  center  of  the  contact 
by  x,=  x0-  xc  where  xc  is  the  distortion  of  the  sample  sur¬ 
face  relative  to  point  O’.  This  distortion  results  from  the 
shear  force  Fs  =  k®  xti  where  K®  is  the  torsional  stiffness  of 
the  cantilever.  For  a  viscoelastic  material  and  oscillatory 


sample  displacement  xjf)  =  A„exp(  icot),  there  is  a  phase  lag 
in  the  response  of  the  contact  so  that  xc  =  A,exp[i(ot/  +  (3)] 
and  x,  =  X,exp[i(cot  -  a)].  (See  the  phasor  diagram.  Fig.  lb). 
The  experimentally  measured  quantities  are  the  amplitude 
and  phase  of  the  tilt  angle  0  of  the  tip,  0  =  x,  IH  = 
0„cxp[  i(  (Ot  -  a)]  where  0„  =  XJH.  The  shear  stiffness  of  the 
contact  is  K,  =  d/'i/dv,. 

The  relaxation  time  x  of  the  polymer  surface  associated 
with  its  viscoelastic  response  is  of  interest  because  it  can  be 
related  to  the  effects  of  intermolecular  constraints  on  seg¬ 
mental  relaxation  [15,16].  We  estimate  x  using  the  simple 
mechanical  Voigt  model  [17].  In  this  model,  the  viscoelas¬ 
tic  component  of  the  contact  is  described  by  a  linear  spring 
of  stiffness  K,  in  parallel  with  a  dashpot  with  damping  co¬ 
efficient  rc.  The  cantilever  is  modeled  by  a  linear  spring  of 
stiffness  in  series  with  the  Voigt  model.  The  stiffness  K,. 
is  related  to  the  effective  storage  shear  modulus  G*  [5].  For 
the  case  of  an  elastic  Hertzian  contact  with  no  slip,  = 
8 G*a  where  a  is  the  effective  radius  of  the  contact  [8,18]. 
This  mechanical  model  is  described  by  the  differential 
equation  F  =  K,.  xc  +  rcdxc/dt  =  K@x,  where  xc  and  x,  are  as 
defined  above.  Force,  F  =  F„exp[;(ft)f  -  y)],  is  related  to  the 
quantities  X,  and  a  by  Fa=  K@X,  and  a  =  y.  The  stiffness  of 
the  contact  is  given  in  terms  of  the  measurable  quantities 
shown  in  Fig.  1  by 
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and  the  response  time  x  by 
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and/ is  frequency,  x,  as  determined  from  eqn.(2),  is  a  par¬ 
ticularly  robust  quantity  since  it  depends  only  on  the  meas¬ 
ured  phase  a  and  the  ratio  of  measured  amplitudes  XtIX„, 
and  does  not  require  knowledge  of  the  cantilever  torsional 
stiffness  K®. 


amplitudes  smaller  than  20  mV  (about  0.5  nm  actual  de¬ 
flection).  Possible  deformation  of  the  tip  [7]  was  ignored. 

The  un-cross-linked  PVE  samples  were  made  from  the 
same  96%  1,2-polybutadiene  material  with  134,000  num¬ 
ber-average  molecular  weight  used  in  a  previous  study  of 
the  bulk  mechanical  properties  [15,16].  From  this  study, 
the  measured  relaxation  time  at  295  K  is  x</  ~  60  (is. 
Young's  modulus  was  measured  by  indentation  to  be  in  the 
range  2-4  MPa.  Films  of  un-cross-linked  PVE  were  cast 
onto  glass  slides  from  toluene  solution,  dried  in  air  for 
about  10-15  min,  followed  by  vacuum  drying  for  another 
10-20  min. 

3.  Results  and  Discussion 

Figure  2  compares  typical  shear  response  data  for 
freshly  cleaved  mica  and  PVE  with  an  applied  load  of  about 
2  nN.  For  mica,  which  is  an  elastic  material,  the  response 
amplitude  X,  increased  linearly  (Fig.  2a)  and  a  was  con¬ 
stant  and  equal  to  0°  (Fig.  2b)  until  a  maximum  drive  am¬ 
plitude  X0max  ~  40  mV  (~2  nm).  At  this  point,  X,  suddenly 
decreased  and  a  increased.  This  is  the  behavior  expected 
when  the  contact  begins  to  slide  [9,23,24].  In  the  case  of 
the  much  softer  PVE,  X,  also  increased  linearly  up  to  125 
mV,  (the  largest  value  studied)  but  with  a  smaller  slope  and 
a  phase  lag  a  =  42°.  The  shear  response  (i.e.,  the  slope  in 
Fig.  2a)  of  PVE  was  measured  from  50  Hz  to  1.2  kHz  and 
the  phase  lag  a  (I's  in  Fig.  2c)  varied  with  frequency/ 
This  frequency  dependent  phase  change  is  direct  evidence 
that  the  contact  is  viscoelastic  [5]. 


FIG  1.  (a)  Schematic  diagram  of  the  SFM  shear  measurement,  (b)  Phasor 
diagram  relating  the  amplitudes  and  phases  of  the  drive  signal  X0,  tip  motion 
Xt,  and  contact  motion  Xc. 


2.  Experimental 

The  SFM  we  used  has  been  described  previously  [19] 
and  is  based  on  the  optical  lever  method.  The  sample  was 
oscillated  laterally  with  respect  to  the  cantilever  by  apply¬ 
ing  a  small  sinusoidal  modulation  to  the  piezoelectric  scan 
tube,  and  the  response  of  the  scan  tube  was  determined 
using  the  double  piezoelectric  response  method  of  Chen 
[20].  Torsional  response  of  the  SFM  cantilever  was  ob¬ 
tained  from  the  photodiode  output  with  a  lock-in  amplifier. 
This  signal  can  be  converted  into  a  shear  force  acting  on 
the  tip  if  the  torsional  force  constant  K®,  and  tip  height,  H 

are  known.  For  the  cantilevers  we  used  [21],  K®  =  (66.6 
N/m)  [22],  and  H  ~  4  pm  as  determined  with  a  scanning 
electron  microscope.  We  use  these  values  in  calculations. 
Xa  was  calibrated  using  measurements  on  diamond,  whose 
shear  modulus  is  about  105  times  larger  than 
poly(vinylethylene),  PVE.  This  avoids  the  need  to  directly 
measure  the  lateral  response  of  the  scan  tube  for  driving 


Drive  Amplitude,  Xo  (mV)  Frequency  (Hz) 


FIG.  2.  Shear  response  (a)  amplitude  X,  and  (b)  phase  a  as  a  function  of 
drive  amplitude  X„  for  mica  and  PVE;  Fn  ~  2  nN.  (c)  Frequency  dependence 
of  phase  lag  a  for  PVE  measured  at  constant  load  (B’s)  and  at  maximum 
load  during  a  force-distance  curve  measurement  (T’s)  such  as  that  shown  in 
Fig.  3. 

The  linear  behavior  of  mica  and  PVE  in  Fig.  2a  is 
somewhat  surprising  from  the  viewpoint  of  contact  me¬ 
chanics.  Micro-slip  near  the  edge  of  the  contact  is  predicted 
to  reduce  Fs  for  a  given  xc  [8,18,25]  and  this  behavior  has 
been  verified  for  macroscopic  contacts  [26].  Peeling  may 
also  be  possible  at  small  displacements  and  would  also  pro¬ 
duce  a  nonlinear  response  [25,27].  Finear  response  is  only 
expected  in  the  limit  of  small  shear  displacements  or  in  the 
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idealized  case  of  a  Hertzian  contact  with  static  friction;  i.e., 
no  slip  anywhere  inside  the  contact  [25].  In  this  case,  Fs  = 
(8 G*a)xc.  The  linear  behavior  shown  in  Fig.  2a  indicates 
that  studies  of  small  contacts  under  shear  should  prove 
valuable  for  testing  the  limits  of  continuum  contact  me¬ 
chanics.  The  contact  radius  for  FN  can  be  estimated  from  a 
Sneddon  analysis  of  the  contact  and  shows  that  the  contact 
radius  is  approximately  equal  to  the  nominal  tip  radius. 

PVE  is  so  compliant  that,  even  for  the  small  forces 
used  here,  the  usual  approximation  that  the  diameter  of  the 
contact  is  much  smaller  than  the  probe  is  not  satisfied. 
Appropriate  expressions  for  an  elastic  indentation  by  a 
spherical  punch  are  given  by  Sneddon  [28]  and  Maugis 

[29] .  We  know  of  no  treatment  for  the  case  of  shear.  The 
contact  radius  for  FN  can  be  estimated  from  a  Sneddon 
analysis  of  the  contact  and  shows  that  the  contact  radius  is 
approximately  equal  to  the  nominal  tip  radius. 

Another  way  to  measure  the  viscoelastic  shear  response 
is  to  apply  the  shear  during  a  force-distance  curve  meas¬ 
urement;  i.e.,  during  the  formation,  loading,  unloading,  and 
rupture  of  a  contact  [4].  Figure  3  shows  data  for  a  freshly 
prepared  PVE  substrate  obtained  at  /  =  1.0  kHz.  The  origin 
for  the  time  axis  is  taken  at  the  point  of  maximum  applied 
load  (C).  The  lower  panel  displays  the  force-distance  curve 

[30] ,  The  sample  was  moved  toward  the  tip  at  constant 
speed  starting  out  of  contact.  Jump-to-contact  occurred  at 
point  A  (t  =  -36  s).  The  approach  continued  until  a  prede¬ 
termined  maximum  load  at  C  (t  =  0  s).  Then  the  direction 
of  motion  was  reversed  until  the  tip  and  sample  separated 
at  the  pull-off  point  E  (  t  =  164  s).  The  rounding  just  prior 
to  pull-off  is  typical  of  both  macroscopic  and  SFM-scale 
viscoelastic  contacts  [31,32];  the  magnitude  of  the  pull-off 
force  and  the  degree  of  rounding  depend  on  the  speed, 
again  indicative  of  the  viscoelastic  nature  of  the  contact. 
For  times  between  points  B  and  D,  the  contact  was  under 
compressive  load;  otherwise,  the  load  was  tensile.  The 
force -distance  curve  was  identical  with  and  without  shear 
modulation  indicating  that  the  applied  shear  modulation 
does  not  significantly  modify  the  development  of  the  con¬ 
tact. 


FIG.  3.  Shear  response  of  a  freshly  prepared  PVE  sample  at  1.0  kHz. 
Lower  panel:  force-distance  curve.  Upper  panel:  amplitude  X,  and  phase 
a  of  tip  motion  on  the  surface. 


The  upper  panel  shows  the  displacement  X,  and  phase 
a  of  the  SFM  tip.  At  jump-to-contact  (A),  X,  jumped  sud¬ 
denly  to  about  1.3  pm,  then  increased  more  slowly  to  about 
2.4  pm  at  maximum  load  (C).  X,  continued  to  increase 
slightly  until  t  ~  60s  even  after  maximum  load  was  reached, 
then  gradually  decreased.  This  decrease  became  steeper 
once  the  force-distance  curve  began  to  round  significantly 
just  prior  to  pull-off  ( E ).  Phase  response  mirrored  the  am¬ 
plitude  response:  a  increased  rapidly  to  about  60°,  then 
decreased  slightly,  reaching  a  minimum  of  about  43°  at  the 
same  time  X,  reached  a  maximum.  For  both  amplitude  and 
phase  response,  the  rate  of  change  during  jump-to-contact 
(at  A)  and  pull-off  (at  E)  is  limited  by  the  overall  response 
of  the  measuring  system.  The  large  decrease  in  phase  just 
after  A  and  the  increase  just  prior  to  E  may  indicate  that 
sliding  occurs  just  after  jump-to-contact  as  well  as  just  prior 
to  pull-off. 

Measurements  on  PVE  were  carried  out  over  the  range 
50  Hz  </<  1.2  kHz  and  were  all  qualitatively  similar  to  the 
results  for  1.0  kHz  (Fig.  3)  except  that  X,  and  a  varied  with 
frequency  due  to  the  viscoelastic  character  of  the  contact 
[5].  Results  were  very  reproducible  for  constant  speed  and 
maximum  load.  The  maximum  X,  and  corresponding 
minimum  a  always  occurred  after  the  contact  was  under 
tension,  even  for  the  slowest  speeds  (50  pm/s).  This  de¬ 
layed  maximum  seems  to  be  a  characteristic  of  viscoelastic 
materials  and  was  not  observed  for  any  of  the  elastic  mate¬ 
rials  we  have  also  studied;  e.  g.,  diamond,  mica,  silicone 
[33],  For  the  data  in  Fig.  3,  X,/X„  «  1,  so  that  eqn.  ( 1 )  can 
be  written  as  a  s  (K0/8G*)(V,/V„)cosa.  Thus,  the  torsional 
response  of  the  force  sensor  is  a  measure  of  the  contact  ra¬ 
dius  a  as  has  been  pointed  out  previously  for  elastic  materi¬ 
als  [6].  The  delayed  maximum  indicates  that  a  reaches  its 
maximum  value  some  time  after  the  maximum  load  has 
been  applied.  Ting  [34]  has  predicted  this  behavior  for 
Hertzian  contacts  in  the  absence  of  adhesion  and  dispersion 
forces. 

The  frequency  dependence  of  a  measured  at  maximum 
load  (e.  g.,  Fig.  3,  point  C)  is  plotted  in  Fig.  2c  (  ’ s). 
These  values  are  slightly  higher  than  the  values  determined 
from  the  linear  response  data  (■'s),  partly  because  the 


Frequency  (Hz) 


Frequency  (Hz) 


Fig.  4.  Voigt  model  analysis  of  (a)  the  response  time  t  and  (b)  stiff¬ 
ness  K  of  the  contact  showing  changes  as  a  function  of  sample  age  for 
PVE. 
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contact  area  had  not  reached  its  steady  state  value  and 
partly  because  the  measurements  were  made  on  a  sample 
that  had  not  aged  as  long,  which  we  discuss  below.  Part  of 
the  difference  may  also  be  due  to  intrinsic  heterogeneities 
on  the  PVE  surface. 

Equations  ( 1 )  and  (2)  were  used  to  calculate  x  and  Kc. 
The  results  are  plotted  in  Fig.  4a  and  b,  respectively,  for 
PVE  as  a  function  of  sample  age,  again  using  data  from  the 
maximum  load  point  in  Fig.  3.  For  the  as-prepared  sample 
(age  <  1  hr),  x  was  about  2  ms  at  50  Hz  and  decreased  to 
about  0.3  ms  at  1  kHz  (Fig.  4,  ’s).  This  reduction  in  x  was 

accompanied  by  an  increase  in  the  stiffness  Kc  of  the  con¬ 
tact:  i.e.,  as  /  increased,  the  contact  became  more  glassy, 
just  as  observed  for  the  bulk  [5].  As  the  sample  age  in¬ 
creases,  x  approaches  the  60  pis  value  observed  for  bulk 
PVE  [16],  After  167  hi'  (■’s),  the  zero  frequency  stiffness 
corresponds  to  G*  ~  2  MPa  and  is  the  same  as  the  bulk 
value. 

4.  Summary 

In  summary,  we  have  demonstrated  the  effects  of  vis¬ 
coelastic  response  to  shear  deformation  in  nanometer-scale 
contacts.  In  particular,  we  show  in  equation  (2)  that 
knowledge  of  the  torsional  stiffness  Kq  of  the  SFM  cantile¬ 
ver  is  not  required  to  determine  the  relaxation  time  x.  As 
the  age  of  the  poly(vinylethylene)  sample  increased,  the 
contact  became  more  elastic  and  the  effective  modulus  and 
relaxation  time  approached  bulk  values.  As  predicted  by 
Ting  for  Hertzian  contacts,  the  maximum  contact  area  oc¬ 
curred  during  the  loading  cycle  well  after  the  maximum 
load  is  applied.  Our  experiments  also  show  that  shear 
modulation  appears  to  provide  an  excellent  way  study  the 
formation  and  breaking  of  small  contacts,  as  well  as  to  test 
the  predictions  of  contact  mechanics  theories  in  the  limit  of 
small  contacts. 
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